Sensitivity of three-body decays to the reactions mechanism and the initial structure 

by example of 6 Be 
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The influence of the initial-state structure and the reaction mechanism on three-body decays is 
investigated by example of the 6 Be continuum populated in neutron-knockout reactions on 7 Be. 
The sensitivity of the 6 Be excitation spectrum and the three-body correlations to the different 
components of the model is demonstrated. It is shown that the spin composition of the initial state 
may have an overwhelming effect on the three-body continuum. The characteristics and structure 
of the second 0^ and 2^ states in 6 Be are predicted for the first time and the conditions for their 
reliable observations are formulated. The effects of interference and the alignment of three-body 
states on the three-body correlations are demonstrated. 

PACS numbers: 25.10.+S, 23.50.+Z, 21.60.Gx, 24.70.+S 
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Introduction. - The study of systems beyond the 
nuclear driplines is an important trend of modern low- 
energy nuclear research. Many of these systems belong 
to the three-body or even few-body continuum and the 
reaction theory for populating these states is not well de- 
veloped. Modern high-precision experiments with exotic 
beams require complicated analyses, advanced theoreti- 
cal treatment, and call for deeper insights in this field. 

Observables in reactions producing unbound systems 
depend on three major ingredients: (i) the structure of 
initial nuclei, (ii) the reaction mechanism, and (iii) the 
final-state interaction (FSI). For very narrow resonances 
(extremely long-lived states), the aspects (i) and (ii) lose 
importance as the structure formed in the reaction has 
enough time to "forget" how it was created. Then for 
a consistent description of the system, it is sufficient to 
study only the decay process (FSI effects) by itself. How- 
ever, exactly when this approach becomes valid is not 
always clear in advance. Clarity in this issue is espe- 
cially important for systems beyond the driplines, where 
the resonant states (often already the ground states) are 
quite broad. 

In this work, we demonstrate the importance of the 
reaction mechanism and the initial-state structure for in- 
vestigations of few-body systems beyond the driplines. 
This is made by example of 6 Be (three-body a+p+p con- 
tinuum) populated in neutron-knockout reactions from 
7 Be projectiles. The first result of these studies has been 
published in Ref. [l[ elucidating the mechanism of demo- 
cratic decay. Good agreement with experimental data 
was demonstrated for both the excitation spectrum and 
the three-body energy-angular correlations over a broad 
range of excitation energy. However, in this compact ex- 
perimental work, many important theoretical issues of 
broader interest by themselves, were left aside. In the 
present work, we would like to focus on the most inter- 
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FIG. 1. Schematic illustrating the 6 Be population formed by 
knocking a neutron out of a 7 Be projectile. 



estingtheoretical results arising from our studies of this 
data [l| . This is a timely message considering the recent 
high interest in the 6 Be system We also think that 

these results should have an important impact on our 
understanding of three-body decays in general and aid 
in formulating experimental strategies for studies of this 
phenomenon. 

Theoretical model. — The major features of the theo- 
retical model were described in Ref. [l[ , but some details 
are needed in the context of this work. The three-body 
a+p+p final-state interactions forming the continuum of 
6 Be are described by solving the inhomogeneous three- 
body Schrodingier equation 
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where ^E^ 4 ") is a wave function (WF) with pure outgo- 
ing asymptotics obtained with the approximate bound- 
ary conditions of the three-body Coulomb problem 0, H| . 
The energy Et is relative to the a+p+p threshold. The 
source term <E> q depends on the vector q of the trans- 
ferred momentum and contains information about the 
precursors and the reaction mechanism. The knockout 
of a neutron from 7 Be is described as a sudden neutron 
removal. After acting with the neutron annihilation op- 
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erator on the 7 Be WF in coordinate space, we obtain 

^ = 1 d 3 r n e^(^ i}ic \^ Bc ). (2) 

Here vector r„ points to the removed neutron, see Fig. [T] 
The 7 Be WF is constructed as an "inert" a-core 
plus p-wave neutron and two protons with coupling 
[(?(^)[(j( 7r i)^( 7r 2)]j]j7Bc- The overlap integral with the 
a particle for this WF is 
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After neutron removal, the terms with coefficients {a, /3} 
lead to the population of the + state in 6 Be while the 
terms with coefficients {7, 6} lead to the population of 
the 2 + state. Within the + and 2 + configurations, the 
ratios a/ j3 and j/S define the spin composition of $ q , 
namely, the probability Ws of configurations with def- 
inite total spin S (coupled spins of the two "valence" 
protons). The components of the source function with 
definite total angular momentum J can be written as 



J,q 



l-r, 2 j\J,S = 0)+r]j\J,S=l), (4) 



where the coefficients r/j, controlling the Ws=o/Ws=i 
ratio, can be expressed via {a,/3,7, 6} using coefficients 
from Table |U 

For single-particle motion, we use the harmonic- 
oscillator WFs whose radial behavior is 

4>{r) = (8/3^F) 1/2 (r 2 Ao /2 ) exp[-r 2 /(2^)]. (5) 

The value ro — 1.8 fm was used in the calculations which 
corresponds within our model to the experimental matter 
radius of 2.31(2) fm for 7 Be. 

The differential cross section is expressed via the flux 
induced by the WF ^( + ) on the remote surface S 
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For the reaction considered, this can be rewritten in 
terms of the density-matrix formalism 



do 



dq dEx dQs 
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TABLE I. The recoupling coefficients from shell-model-like 
(jj) coupling in the source to (Is) coupling of the three-body 
model. 

j 1 ,32,J: 3/2,3/2,0 1/2,1/2,0 3/2,3/2,2 3/2,1/2,2 
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FIG. 2. Excitation spectra of (a) + and (b) 2 + (b) states 
as function of the spin composition of the source function. 
The magenta-colored marked grids correspond to the first and 
second excitations and , see Fig. [4] Arrows show the 
positions corresponding to definite shell-model configurations 
in the source. 



Some ingredients of the Eqs. (|6I7I) are illustrated in Fig.[T] 
In Eq. ([7]), the contributions of the three-body dynamics 
(amplitudes Ajm) and the reaction mechanism (density 
matrix Pj^f ) are explicitly separated. For direct reac- 
tions, the density matrix has an especially simple form in 
the frame with the z axis coinciding with the direction of 
the transferred momentum q. We use two limiting forms 
of the density matrix: 
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P2M — 1/5 j P20 



Poo 



Poo — 1 ' P20 — 1 ■> 



P2I = POO = COS(0 2O 



(<M/\/5, (8) 
(9) 



In the sudden- removal model of Eq. ([2]), there is no align- 
ment of the final state. This should lead to Eq. (jHJ with 
the relative phase of the + and 2 + states (^20 = 0. How- 
ever, it is clear that some alignment should be introduced 
by a realistic reaction mechanism. To check our sensitiv- 
ity to this, we also used the density matrix of Eq. ([9]) cor- 
responding to the completely aligned case and also kept 
020 as a parameter in the both above cases. In order 
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FIG. 3. Left axis: theoretical 6 Be excitation spectrum to- 
gether with the results of the MC simulations superimposed 
on the data Right axis: predicted evolution of the spin 
composition in 6 Be for + , 2 + , and for the total spectrum 
fitted to experimental data (the sum of + and 2 + contribu- 
tions). 



to compare with experiment, Eq. ([7]) provides theoreti- 
cal input for Monte Carlo (MC) simulations used to deal 
with the bias introduced by experimental apparatus. 

The reaction model and the initial-state WF can be 
seen as simplistic. However, this fits the aim of this pa- 
per to provide a "proof-of-concept" demonstration, where 
interdependencies on different aspects of the model are 
as transparent as possible. 

Sensitivity of the 6 Be spectrum to Ws=o/Ws=i is 
demonstrated in Fig. [5] For some range of parameters, 
the conventional picture of 6 Be excitation is observed 
with the 0+ and 2+ states at 1.37 MeV and 3.05 MeV, re- 
spectively, with significant variations taking place only in 
the high-energy "tail" of the distribution. However, with 
other variations of the parameters, the "normal" states 
of 6 Be "fade" and even completely disappear, while new 
broad peaks arise for these two states at ~ 12 MeV and 
~ 7 MeV, respectively. Thus, our reaction model links 
variations of the 6 Be excitation spectrum to variations of 
the structure of 7 Be. In this work we would like to investi- 
gate this link quantitatively without reference to realistic 
structure of 7 Be in order to find boundaries for possible 
scale of effects and establish principal opportunity to use 
three-body-continuum spectra as tools to study the spec- 
troscopy of the precursor. 

The broad range of spectra obtained for + and 2 + 
continuum provide an opportunity to fit the experimen- 
tal spectrum of 6 Be. This fit, obtained in [l|, is shown 
in Fig. [3] together with the MC and experimental data. 
It corresponds to the simple case of pure 5 = pop- 
ulation in the source function {rjj = 0) correspond- 
ing to a parameter set for initial state of {a, /3, 7,(5} = 
{0.42,0.3,0.49,0.7}. It should be noted that the spin 
content of the final state may have nothing in common 
with the spin content of the source function as the spin 
quantum number is not conserved by the Hamiltonian 
([T]). The predicted final- state spin composition (Fig. [31 
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FIG. 4. (Color online) Excitation spectra and spin compo- 
sition for the 6 Be continuum explicitly representing the first 
and second resonances for (a) = + and (b) 2 + . Note, 
the broken abscissa in (a). 



right axis), evolves rapidly between the location of the 
+ peak and past the location of the 2 + peak indicating 
important modifications of the nuclear structure in this 
energy region. 

Second + and 2 + states — The above observations 
allow us to determine the properties of the 0^" and 2^ 
states. We just need to choose the parameter settings 
minimising population of the normal resonant peaks. An 
important feature of "pure" first and second states, illus- 
trated in Fig. H] (see also arrows in Fig. [2]), is that the spin 
composition depends weakly on energy. The existence of 
a common spin structure independent on energy, allows 
us to interpret and as different states, although 
they can be represented by (relatively) broad overlapping 
structures. This property of "pure" states is in sharp con- 
trast with the spin evolution for the "composite" situa- 
tion of Fig. [3] The spin-content ratios 0]*~/0^ and 2+/2^~ 
provide a simple structural idea: these states are part- 
ners, in the sense they are the orthogonal combinations 
of the 5* = and 5 = 1 configurations. E.g. for 2 + states 
it can be seen that the Ws=o/Ws=i ratio is 2:1 for the 
first and ^1:2 for the second resonance. It can be seen 
from Figs. [5] and S] that the "pure" states are best popu- 
lated from source functions which are close to pure shell 
configurations in Eq. (PJ and in fact their structure is 
reasonably close to such pure shell configurations. 

The above insight on the structure of the 0^ and 2^~ 
states of 6 Be provides clear guidelines for experimental 
searches: the reaction mechanism should enrich the S = 
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as being completely denned by the dynamics of the three- 
body motion itself. The parameters 
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FIG. 5. (Color online) Dependence of the 6 Be excitation spec- 
tra on the radius parameter ro of the 7 Be WF. All spectra are 
normalized to a maximum value of unity. 



1 component in the final state. 

Radial dependence. - The sensitivity to the radial 
characteristics of the source is presented in the Fig. [5] 
It is practically nonexistent for 0^ : variations take place 
only in the "tail", a few decay widths higher than the 
resonance position. The l\ sensitivity is quite small: the 
resonance width is affected on the level of 15%. However, 
with increasing excitation energy of the states (and hence 
with increasing decay width) the effect grows. There is 
about 0.8 MeV uncertainty of the 2^ position connected 
with the radial extent of the source, and the profile of 
the cross section is strongly affected. Variation of 0^ 
properties is so large (few MeV in peak position) that 
that the properties of such a state cannot be discussed 
without a detailed account of the reaction mechanism. 

Three-body correlations. — The three final momentum 
vectors in Eq. (J6)) and Fig. [1] can be transformed to the 
transferred momentum q and the Jacobi momenta k x , 
k y in Eq. ©: 



A 2 ki - Aik 2 



A 3 (ki 



(Ai + A 2 )k 3 



A x + A 2 + A 3 
E T = E X + E V = k 2 x /2M x + k 2 y /2M y , (10) 

where M x and M y are the reduced masses of the X and 
Y subsystems, see Ref. @ for details. 

The orientation of q is not a dynamical variable of the 
model Eq. ([2]) and thus is not present in Eq. ([7]). The five- 
dimensional "hyperspherical solid angle" £1$ includes two 
degrees of freedom describing the "internal correlations" 
of the three-body system which are ordinarily considered 
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h) = (k a • ky)/(k x ky) , (11) 



provide what we call "complete energy-angular correla- 
tions" . For &3 — > k a we get the "T" Jacobi system, where 
e describes the energy correlation in the p-p channel. For 
&3 k Pi , the correlations are obtained in one of two pos- 
sible (i = 1, 2) "Y" Jacobi systems, where e describes the 
energy correlation in the a-p channel. 

The other three degrees of freedom (Euler angles) de- 
fine "external correlations" as they describe the orienta- 
tion of the three-body system as a whole. Correlations 
for the "external" degrees of freedom are evidently de- 
fined by the reaction mechanism. 

Alignment and interference. - The internal three- 
body correlations for excited states cannot be separated 
from external information in any practical experiment. 
The excited states typically have nonzero J making align- 
ment possible and they reside on the "tails" of the lower- 
energy excitations so that their amplitudes can interfere. 
From a theoretical point of view, the inclusive energy 
£ distributions should be free of interference effects and 
experience tell us that the angular cos(0^) distributions 
are only weakly affected. However in experiments, the 
bias of the apparatus introduces cut-offs and distortions 
which may induce correlations via loss of orthogonality 
for configurations with different angular momenta. Such 
"induced correlations" are specific for the experiment and 
should be evaluated by careful MC studies. 

A demonstration of alignment/interference effects on 
experimental results is provided in Fig. [6] As it is clear 
that the effects are largest for strong + /2 + mixing, we 
have chosen the Et = 2 — 3 MeV energy range of the data 
1] . MC results for the two limiting cases of completely- 
aligned ((9l) and "isotropic" (HJ) density matrices, each 
with three different phase settings 020 = {0,7r/2,7r}, are 
shown. In Ref. [l[ , we used the settings of Eq. © with 
020 = which is consistent with the data and corre- 
sponds to the model of Eq. @. However, at the mo- 
ment we cannot exclude that reality is different: further 
analysis of the data is needed. Fig. [5] shows the energy 
distribution in the Jacobi "T" system and the angular 
distribution in the Jacobi "Y" system which appear to 
be the most sensitive of the correlation observables. The 
scale of local variations in the MC distributions is about 
20%. The curves are well separated and some deviate 
considerably from the experimental results. 

It should be understood that all the MC curves cor- 
respond to the same theoretical distribution and their 
variations are due to the bias introduced by experimen- 
tal setup. Thus we conclude that in analysing high- 
precision correlation data for excited states where inter- 
ference/alignment effects become possible, a consistent 
treatment of the reaction mechanism becomes inevitable. 
The effect depends strongly on the quality of the experi- 
mental setup; it should vanish for an "ideal" instrument. 
The moderate level of variations in Fig. [5] is connected 
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FIG. 6. (a,b) "Internal" energy distributions in the Jacobi 
"T" system and (c,d) angular distributions in the Jacobi "Y" 
system. Upper and lower rows of panels correspond to dif- 
ferent Et ranges. The experimental data are shown by the 
hollow diamonds. Theoretical curve (thick, solid, red) is given 
in each panel with some offset to simplify a perception. MC 
curves for different alignment /interference settings are ex- 
plained in the legend. 



with the very high efficiency of the setup in Ref . [l[ . On 
lower-quality setups, alignment /interference effects can 
produce very large and poorly controlled modifications. 
This could be a part of the explanation for the strong 
deviation of the data of Ref. Q from the other recent 
experimental studies 0, d, @ . 

On existence of higher- energy negative parity excita- 
tions. - The availability of higher excitations in the 
spectrum of 6 Be [l[ may drastically affect the proposed 
interpretation of the data. In Ref. [6j , a strong population 
of presumably {0~, l - , 2~} states in the charge-exchange 
reaction induced by AL = 1 transitions was observed 
above the first 2 + state. We were cautious about the 
population of such negative-parity states in the data of 
Ref. [![. A dedicated search was performed for asym- 
metries in certain distributions due to the interference 
of positive and negative states. However, no significant 
indication of these asymmetries was found. Also, the 
dominating shell structure of the 7 Be precursor does not 
imply a strong population of negative-parity states in 6 Be 
within the sudden-removal model. Thus we have confined 
our interpretation of the 6 Be continuum with just the + 
and 2 + states. 

Previous theoretical results on 6 Be. — Three-cluster 
calculations (microscopic and three-body) for 6 Be have 
been performed in a number of studies 0412 ■ In Ref- [3 , 
the energies and widths {Er,T} of the 0^" and l\ res- 
onances were predicted to be {3.5,6.1} and {5.2,5.6} 
MeV. These values are quite different from our predic- 
tions {~ 12,- 14} and {~ 7,~ 7} MeV. The reason for 
this discrepancy is easily understood: our prediction is 
that the second and 2^~ states are "spin-complements" 
of the first 0^ and 2+ states. States of such a nature could 
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FIG. 7. Comparison between the energy distributions for the 
2 + resonance of 6 Be in "T" and "Y" Jacobi systems obtained 
in the present work and from Ref. 



not have been obtained in the calculations of Ref. [9] as 
the S = 1 component was omitted in their model space. 

Among these other theoretical studies, only in Ref. [l2[ 
are theoretical correlations (energy distributions for 2 + ) 
presented. The predicted correlations for the + and 2+ 
states (presumably related to [13]) are provided in 
but in a form which make them difficult to interpret as 
only MC results specific for experimental setup of that 
work are shown. The distributions of Ref. are not in 
agreement with our results (see Fig. [7]) and hence with the 
recent highly-accurate data of Refs7[l], 0] . This observa- 
tion sheds doubts on applicability of the methods of [l2[ 
to Coulombic three-body decays in general. Our work 
imposes a new standard of a sophistication required from 
theoretical calculations to analyse modern high-precision 
data including correlations. 

Conclusions. — This work provides important quali- 
tative insights into the question of which aspects of three- 
body decays can be understood based on the dynamical 
description of the final state alone and which also require 
an adequate treatment of the initial state and the reac- 
tion mechanism. The major and novel observations of 
this work are the following. 

(i) The population of the 6 Be Of ground state is very 
stable to variations of the initial-state structure and the 
reaction mechanism. However, the latter effects become 
increasingly important with increasing excitation. No 
sensible description of the continuum above 5 — 7 MeV 
can be given without their proper treatment. 

(ii) The excitation spectrum of the three-body contin- 
uum of 6 Be up to 15 — 20 MeV is found to be very sen- 
sitive to the spin composition of the source function in 
Eq. (JTJ) . In our calculations, this was parametrized in 
terms of the spin content of precursor. Thus we find that 
three-body decay can be used as a sensitive instrument 
of nuclear spectroscopy if the reaction mechanism is well 
established. 

(iii) A procedure to identify the second 0^ and l\ states 
of 6 Be is proposed. These states are found to be "spin- 
complements" of the well-known Of and 2^ states, pro- 
viding a guideline for their experimental observation. 

(iv) Alignment and interference effects are observable in 
experimental data due to the unavoidable experimental 
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bias. They were found to have an important impact on 
the measured three-body correlations. Proton-proton en- 
ergy distributions are especially affected. This indicates 
that caution is needed in studies of N-N correlations in 
decays of excited (or/and broad) three-body states in 
general. 
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